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Abstract 

In this paper we analyze a constant multiple-input multiple-output interference channel where a set 
of active users are cooperating through interference alignment while a set of secondary users desire 
access to the channel. We derive the minimum number of secondary transmit antennas required so 
that a secondary user can use the channel without affecting the sum rate of the active users, under 
a zero-forcing equalization assumption. When the secondary users have enough antennas, we derive 
several secondary user precoders that approximately maximize the secondary users' sum rate without 
changing the sum rate of the active users. When the secondary users do not have enough antennas, we 
perform numerical optimization to find secondary user precoders that cause minimum degradation to 
the sum rate of the active users. Through simulations, we confirm that i) with enough antennas at the 
secondary users, gains equivalent to the case of all the users cooperating through interference alignment 
is obtainable, and ii) when the secondary users do not have enough antennas, large rate losses at the 
active users can be avoided. 



I. Introduction 

Interference alignment (IA) [|2), is an interference management technique broadly described 
as confining the interference to a subspace of the received signal space such that an interference 
free subspace becomes available for the desired signal. Unlike techniques such as ignoring the 
interference |ffl, decoding/cancelling the interference [5] or avoiding the interference (orthogonal 
access), IA achieves the maximum degrees of freedom of the K-user interference channel. 

Interference can be aligned using signal levels ||6), time or frequency channel extensions 
[|2j, and multiple antennas at the interfering nodes (7J-||9j. IA using the signal levels (lattice 
alignment) requires a baseband module not compatible with near future wireless networks |10j. 
IA over time extensions of the channel is a mere theoretical tool requiring non-causal knowledge 
of channel state information. IA over frequency extensions of the channel, in its original form 
0, requires a large number of channel uses to fully exploit the power of IA, and is not suitable 
when global frequency synchronization is not available. Therefore, aligning the interference over 
the spatial dimension, i.e. IA over constant multiple-input multiple-output (MIMO) channels, is 
the most attractive physical layer signaling approach. In this paper, we focus on MIMO IA. 
In general, a high computational [11|, [ fT2| and overhead [13|, [14| cost is associated with 



finding the IA precoders and combiners. Therefore, once a set of users have aligned their 
interference, it is desirable to retain their alignment status until the channels change. Most future 
interference-limited wireless networks, however, will be packet-switched with bursty data traffic, 
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requiring frequent changes in the number of active users JTOj. Also, the feedback overhead p"5| 
and the number of antennas at each node Q practically limit the number of user pairs that can 
cooperate through IA. Therefore, in an IA network, there will be nodes that are arriving to the 
network and, although not included in the existing IA setup, wish to communicate with their 
receivers. Determining when such nodes can be admitted to the network (are allowed to transmit) 
and/or developing transmission strategies for them is referred to as user arrival. Although there 
exist a few prior strategies for user arrival and precoder design for MIMO networks, they either 
do not consider IA or they are restricted to opportunistic access methods p6j-^9j. 

In this paper, we consider a set of active users exhausting the network resources by coopera- 
tively utilizing IA and a further set of secondary users who wish to communicate in this network. 
We assume the secondary users are required to have minimum impact on the performance of the 
active users, defined as the sum rate of the active users after a zero-forcing (ZF) receiver, and the 
active users ignore the presence of the secondary users when designing their precoders and/or 
equalizers. We first compute a zero-impact threshold for the number of secondary transmitter 
antennas where secondary transmitters with more (or equal) antennas than this threshold can use 
the communications medium without degrading the sum rate of the active users. Then we find 
optimum and suboptimum secondary user precoders for two cases: (i) when the secondary users 
satisfy the zero-impact threshold and (ii) when they do not. 

When the secondary users satisfy the zero-impact threshold, there exists a set of precoders 
that will not degrade the sum rate of the active users. Thus, the secondary users can optimize 
an objective function of their own link by selecting a precoder from this set. We choose the 
achievable sum rate as the objective function and derive optimum and suboptimum precoders for 
the case of one and two secondary users respectively. For more than two secondary users, we 
choose the degrees of freedom (DOF) as the objective function, defined as the slope of average 
sum rate (b/s/Hz) versus logarithm of signal-to-noise ratio (dB) at high transmit power. We 
propose successive IA precoding and show it is optimum for various network setups determined 
by the number of active users, secondary users, and antennas at each node. When the secondary 
users have fewer antennas than the zero-impact threshold, we search for secondary user precoders 
causing minimum degradation to the sum rate of the active users through a steepest descent 
search over the Grassmann manifold. For this numerical optimization, we propose three initial 
solutions of varying degrees of complexity. Our initial work in [ 1 1 only deals with the case of 
enough antennas at the secondary users. This paper elaborates on the claims of [[TJ, provides 
new analysis for the case of enough antennas at the secondary users, considers the case of not 
enough antennas at the secondary users, and presents new simulations results. 

The remainder of the paper is organized as follows. In Section [TT] we present the system 
model and review IA in the f^-user constant MIMO channel. Admission of the arriving users 
is considered in Section IHlJ Section [TV] presents precoder designs for the secondary users with 
large number of transmit antennas. In Section [V] we find secondary user precoders when the 
secondary users do not have enough transmit antennas. Numerical results are presented in Section 



VI followed by concluding remarks in Section VII 



Notation: Upper and lowercase bold letters represent matrices and vectors. A*, A T and A -1 
are conjugate transpose, transpose and inverse of A. A(n,m) is the element on the nth row 
and the mth column of A. tr(A) and rank(A) are trace and rank of A. || • \ \ F and || • || are 
the Frobenius and zero norms. [A,B] is the matrix constructed by horizontal concatenation of 
matrices A and B. \a] is the smallest integer larger than or equal to a. Kronecker product is 
represented with cg>. For two matrices A and B with dimensions a\ x a-i and b\ x b 2 , A © B 
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is a block-diagonal matrix of dimension (a x + &i) x (a 2 + b 2 ) with matrices A and B on its 
main diagonal. Ijv is the N x N identity matrix, aX b is the matrix of dimension a x b with all 
elements equal to zero. Finally, O is the set of all matrices with orthonormal columns. 



II. System Model 

Consider a -fd-user MIMO interference channel shown in Fig. [T] where the ith transmitter 
and receiver pair are equipped with Mj and iVj antennas. Each transmitter % uses a precoding 
matrix Fj of dimension M\ x di to transmit di streams to its corresponding receiver. In this 
paper we assume Fj 6 O. Note that precoders with orthonormal columns are useful for limited 
feedback codebook design [20], in many cases unitary precoders are optimal when there exist 
a peak power constraint [|21[, and waterfilling over the eigenmodes of the effective channels 
can always be used to further improve the performance of the system after appropriate unitary 
precoders have been found. For each time instant, the received signal at the ith receiver, with 
perfect timing and synchronization, is 

K a 

Yi = ^ H ifc FfeX fe + Zi i = l,...,K (1) 

fc=i 

where H ik is the matrix of channel coefficients of a block fading channel between transmitter k 
and receiver i, the transmitted signal from the ith node is x, with power constraint E{x*Xj} = P 
and Zj is the AWGN with elements in CjV(0, a 2 ) where a 2 is the noise power spectral density 
which includes thermal noise and white excess interference from unaccounted sources. Note that 
similar to j9), we could incorporate colored noise by prewhitening or modifying the proposed 
algorithms but this is not considered in this pap er. 

A fT a -user system of IA is feasible if there exists a set of matrices W = {Wi, . . . , W^ a } 
such that, given the system model of ([T|), the following constraints are met Q 

f rank (W^F*) = di y . , R , 2 

\ W,H ifc F fc = VA; ^ i Vz ' fc G 1 X ' ' ' ' ' A(lj ' {Z) 

where di is the number of interference-free streams the ith active transmitter can send to its 
receiver. The linear equalizer presented in (221 and the projection matrix presented in |9l Section 
III.A] are examples of a possible receive filter Wj in ([2]). We assume Hik has full column 
rank for all i,k and the set of {Mj, iVj, dj, K a }, Vz, is feasible |7j. Through IA, interference 
from the active transmitters at the ith active receiver is confined to an iVj — di dimensional 
subspace. We denote the non-unique basis for the interference subspace as Cj. We assume the 
ith active transmitter sends di streams and each active receiver uses a ZF equalizer, given by 
Wj = 0di,Ni-di] [HjjFj,Cj] _1 p2| . The total achievable sum rate of the active users is |22| 



^ u u m = Y,Ym [ 1 + _n . ] , o) 



"sum / ■ / ■ 
i=l n=l 



((FIHIP.H^F,)- 1 ) e; 



where Pj = (1^ — CjC*) is the projection matrix into the nullspace of the interference subspace 
at the ith active receiver, 7 G = ^ and e n is the nth row of I n . Note that the assumption of 
a ZF receiver at the active receivers is for mathematical tractability and developing intuition 



into the problem. The design of section V-B is independent of the receiver type and the design 



of Section V-A can be extended to any receiver type that satisfies the rank constraint in |2j 



We assume global perfect CSI knowledge and that the secondary users can use the CSI of the 
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active users' network to find the precoders/equalizers used by the active users. Although this 
assumption is not practical, the performance achieved by the developed techniques in this paper 



can be used as a benchmark for the case that CSI is relayed via limited feedback [13|, [14|. 



III. Secondary User Admission 

Assume K s secondary users request access to the interference alignment group. Define K = 
K a + K s , and let the secondary users be the last K s users in the ordered set of user indices 
K, = {1, . . . , K a , K a + 1,...,K}. The received signal at the ith < K a receiver is y, = 
Yl^li Hj fc F fc x fc + Yld=K +1 H^F^x^ + Zj. We assume the first K a users do not change their 
precoders or receiving filters when the secondary users join the network. This has two important 
implications: a) the active transmitters do not help the transmission of the secondary users, b) 
the active receivers are not aware of the interference from the secondary users. 

In the presence of the secondary users, the sum rate of the first K a users changes from Q to 

K a di I ,, \ 

Ki^yy^g i+— , ^ s— , (4) 

where we have assumed that H ik is independent of H nm for i ^ n or k ^ m. Under the 
assumptions that the active users are ignorant of the secondary users, a reasonable objective is 
for the secondary users to not degrade the sum rate of the first K a users. Comparing Q and ([3]) 
implies that ^2 k=Ka+l ^WjH ifc F fc F£H* fc W* should ideally be zero or equivalently 

W,H 4fc F fc = ie{l,...,K a },ke{K a + l,...,K}. (5) 

It can be inferred from ([5]) that the interference from the secondary users should be confined to 
the interference subspace of the active users p2l. Therefore, by requiring that PjHjfcFfc = for 
i G {1, . . . , K a } and k e {K a +1, . . . , K}, the precoder of the kth transmitter should satisfy 

H fe F fc = (E Xc, iJV . )xdfc , (6) 

where = [(PiHi^)* , . . . , (P^H^)*]*. Using ([6]), Lemma [T] provides the required minimum 
number of antennas at each secondary user, called the zero-impact threshold, which enables the 
secondary nodes to communicate, potentially, without affecting the sum rate of the active users. 

Lemma 1. The kth secondary user can transmit d k streams without degrading (|3]) if 

K a 

M k >J2di + d k k G {K a + 1, . . . , K}. (7) 

i=l 

Proof: As H ik has full column rank and Pj is a projection matrix into a di dimensional 
subspace (rank(Pj) = di), rank(PjH ifc ) = di p3| Property 3.9.b]. Hence, an F fc ^ satisfying 

K a 

^ can be found if M k > rank(LL) = £ di. As each secondary transmitter sends d k streams 

to its receiver, F k (the nullspace of H fc ) should have rank of d k and the result follows. ■ 
Based on Lemma [TJ we divide the technical contributions of this paper into two parts, 
i) The number of transmit antennas at each secondary user satisfies Q. In this case, the 
secondary users can also optimize some performance criterion for their own link. 
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ii) The number of transmit antennas at each secondary user does not satisfy ([7]). The secondary 
users once admitted degrade the sum rate of the active users. 
In the next two sections, we derive optimum and suboptimum secondary user precoders for each 
scenario based on K s and the available transmit antennas at the secondary users. From this point 
on, for mathematical tractability and to simplify the exposition, we assume the first M a x iV a 
K a users each transmit d a streams and the last M s x N s K s users each transmit d s streams. 

IV. Secondary Users Not Interfering with The Active Users 

Based on (|6]) and Q, the secondary transmitters need at least K a d a + d s antennas in order not 
to decrease the sum rate of the active receivers. This means each secondary user has to generate 
a precoder of rank d s using M s — K a d a bases. Therefore, when M s — K a d a > d s , the secondary 
users can use the extra degrees of freedom in designing their precoders to not only satisfy ([6]) 
but also improve the performance of their own link. For the rest of this section, we assume 
M s > K a d a + d s and that the secondary users want to maximize their sum rate given by 

Kim = E lo S det f 1 + jK kk F k FlW kk I-A , (8) 

k=K a +l ^ s ' 

K 

where X k = (l]y k + Yl ^HfcjFjF*!!^) captures the interference from all the transmitters. 

Let the right singular vectors of spanning its nullspace be the columns of V fc . It is inferred 
from (|6]) that the columns of F k have to be a linear combination of the columns of V^; i.e., 
Ffc = VfcGfc, for some matrix G k . As both F^ and have orthonormal columns, columns 
of Gfc ^are also orthonormal, i.e. G O. Using ((8]), the optimum secondary user precoders, 
Ffc = V fc Gfc, are found by solving 



argmax V log det I + ^H fcfc V fc G fe G^V*H* fc X^ , (9) 



k=K a +l 

s.t. H fc V fc G fc = k = K a + 1,...,K. 



Solving ([9]) is difficult in the most general case; we focus on several important special cases. 



A. Single New User 

In this section we assume a single secondary user is arriving to the network. Then the optimum 
secondary user precoder solving (|9]) can be found in closed-form and is given by Lemma [2} 

Lemma 2. For K s = 1, the columns of G k solving ^ are the d s most significant eigenvectors 
ofViW kk I- k l ll kk V k . 

Proof: Let the summands in the objective function of (|9]) be f(G k ). Using the identity 
det (I+AB) = det (I+BA) for any two matrices A and B, rewrite f(G k ) as 

f(G k ) = log det (i + gvjHJ^HttVfcGfcG^ . 

Let the eigenvalue decomposition of \lH.* kk I k 1 ii kk V k be Ujfc£ fc U£ where S fc is a diagonal 
matrix holding its eigenvalues. Then f(G k ) = log det fl+ ^S fc U^GfcG^,Ufc j , and using the 
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for K s > 1 is not 



Hadamard inequality, the G k diagonalizing the term inside the det maximizes f(G k ). A solution 
is found by setting the columns of G k to the d s most significant eigenvectors in U^. ■ 
If the secondary user was not optimizing its own link, G& could be set to a random uni- 
tary matrix. Call this approach random orthonormal design. Moreover, the secondary user 
could ignore the IA nodes and maximize its own achievable rate by adopting a selfish pre- 
ceding through setting the columns of its precoder to the d s most significant eigenvectors of 

H*kk \ In. + Y a S H KiFiF*H* K ^j H KK . Note that generalizing Lemma 
trivial while random orthonormal and selfish precoding can be used for arbitrary K s . 

B. Two New Users 

Now assume K s = 2 where the goal of the secondary users is satisfying ([6]) while maximizing 
the total achievable sum rate of the secondary users given by ([8]). In this case, I k in ([8]) is 

K a 

Xfc = I +?E H ^ F ^ H ^ + Z Hfc " F ^ H ^ <l,k£{K-l,K}, (10) 

where q^k. Corollaries [T] and [2] give approximate solution for ([9]) in some special cases. 
Corollary 1. When K a 3> (K s — 1) > 1, (|5J> w solved by setting the columns of G k to the d s 
most significant eigenvectors of V* k H* kk (l Na + g ESl H fc; F i F * H fci) H fefc V fe . 



Proof: When K a ^> (K s — 1) > 1, the second term in ( flO] ), representing the interference 
from the active transmitters, dominates the third term which is the interference from the other 
secondary user. Ignoring the third term in ( [10] ) changes the optimization of ^ to the case of 



K s — 1 which is solved using Lemma [2] and the result follows. 

In Corollary [T| the secondary users can independently design their precoders. A higher sum 
rate is obtainable, however, by successive application of Lemma [2] to the nodes in the ordered 
list of {K a + 1, . . . , K} such that the kth secondary user ignores the interference from the 2th 
transmitter for i £ {fe + 1,..., K}. Corollary [2] gives a more precise description of this approach 
where for conciseness we ignore the gain attainable through ordering the users. 

Corollary 2. When K a ^> (K s — 1) > 1, Q is approximately solved by computing Gk-i using 
Corollary [7] and setting the columns of Gk to the d s most significant eigenvectors of 

-l 



V*kH*KK ( In s + ^^K{K-l)F K-lF*K--i^*K{K-r) + H ^ F i F i H Xi ] HkK^K- 

In an IA network configured for K a users, the number of secondary users arriving is expected 
to be less than the number of users currently in the network which justifies the assumptions 
of Corollaries [T] and |2j Moreover, both Corollaries [T] and [2] can be readily extended to more 
than K s = 2. As K s increases, the interference from the other secondary user(s) becomes non- 
negligible and we expect the performance difference between Corollaries [T] and [2] to increase. 

In Corollary |2} one of the secondary users interferes with the other one. A more involved 
successive application of Lemma [2] can confine the interference from the Kth. transmitter at the 
(K — l)th receiver to the (iVj^_i — Ylf=i di) dimensional subspace spanned by the interference 

K a K a 

from the active users which implies M K > + {N K -i — di) + d K = N K _i + d K . 

t=i i=i 



K a 
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C. More Than Two New Users 

Solving (|9]) for K s > 2 is equivalent to solving the general sum rate maximizing precoder of 
the MIMO interference channel and, to date, a closed-form solution directly solving it does not 
exist JT6). As an alternative, we provide a precoder design maximizing the pre-log factor of ([8]) 
at asymptotically high transmit power for a certain (and plausible) network configuration. 

From ([6]), the ith secondary transmitter has to allocate K a d a of its spatial dimensions to pre- 
align the interference it causes to the active receivers by selecting its precoder in the column 
space of Vj. Interference from the active transmitters at the kth secondary receiver spans a 
K a d a dimensional subspace given by the interference basis of {HfciFi, . . . , Yl kKa F Ka ). Let 
columns of W fc span the left nullspace of [(H fc iFi) T , . . . , (H kKa F Ka ) T ~\ . An equalizer at the 
A;th secondary receiver post multiplied by always cancels the interference from the active 
transmitters. Let the channels between the secondary users, H k i for k, i e {K a + 1 K}, be 
replaced by (M s — K a d a ) x (N s — K a d a ) dimensional effective channels of WfcHfcjVj. Then the 
networks of the secondary and active users will become disjointed meaning that the secondary 
nodes will not cause/receive any interference to/from the active nodes. 

Assume M s = N s and let M s = N s = M s - K a d a . If K s < 2^ - 1 R, or M s = 
(when d s = 1), we can perform another level of IA between the secondary users over their 
effective channels. Hence both the active and the secondary users will be performing I A among 
themselves and while the interference from the secondary transmitters is pre-aligned to the 
interference subspaces at the active receivers, the interference from the active transmitters is 
always canceled at the secondary receivers. We call this method successive IA and conjecture 
that in special cases, it achieves the same DOF as if all the K users had performed IA together. 

Conjecture 1. Consider a 3K-user interference channel for K G Z + where the transmit- 
ter/receiver pairs are divided into K groups of 3 users each, {G±, G2, ■ ■ ■ , Gk }, such that 
the nodes of the ith group have 3z — 1 antennas. Performing successive IA on the kth group, 
2 < k < K, through creating effective channels between the nodes of G k based on the 
interference subspaces and precoders of the {l,...,k — 1} groups achieves the same DOF 
as if all the 3K users had done IA together. 

Using counting arguments, it is possible to show that successive IA achieves 3K DOF. Proving 
that IA can only achieve 3K DOF, however, is not trivial. We have numerically confirmed (up 
to K = 5) that only 3K DOF is achievable if all the nodes perform IA together (for DOF > K, 
the algorithm of [|9jSection III. A] does not converge). We also use |7| Theorem 2] to further 
support Conjecture [T] Let all the nodes except the Ah node (in the ith group) transmit a single 
stream. Then, the number of variables (N v ) and the number of equations (N e ) in ([2]) are 

3K 



N V = Y^ dk{M s + N k - 24) = 9K 2 - 3K - 61 + 4 + d e (6i - 2 - 2d e ), (11) 

k=l 
3K-1 

N e = ^ (de + 3K-2)+d e {3K - 1) = {3K - 1) (3K - 2 + 2d e ) . (12) 



k=l 



From ( flT] ) and ( |T2| ), since di > 1 and i < K, this system is always improper (N v < N e ) 

N v = 9K 2 - 3K - 6i + 4 + d £ (6i - 2 - 2d t ) < N e = (3K - 1) (3K - 2 + 2d t ) 
6d e (i — K) — 2d] - 6(z + K) < -2. 



8 



Note that if K s > 2^ — 1 or M s ^ [^^] (assuming M s = N s ), successive IA is not possible 
and, in general, the secondary users precoders have to be found through numerical optimization 
methods. Secondary users, however, can use suboptimal precoding designs including: 

• Each secondary user independently designs its precoder based on Corollary [TJ 

• Secondary users successively design their precoders based on Corollary [2j 

• Performing successive IA on a subset of secondary users and using any combination of the 
previous suboptimum methods for the rest of the secondary transmitters. 



We present numerical comparison of these methods in Section VI and leave analytical comparison 
of such suboptimum solutions for future work. 

V. Secondary Users Interfering with the Active Users 

In this section we assume that the number of transmit antennas at the secondary users is less 
than the zero-impact threshold given by ([7]). This implies that (|5]) cannot be satisfied and the 
secondary users are bound to change the sum rate of the active users. Consequently, we desire 
the secondary users to cause minimum degradation to the sum rate of the active users given by 
(Q. Defining Q; = (F* H^PjHjjFj) -1 , we seek a solution to the problem of 

argmax VVlog 1 + — - r r . (13) 



Like the sum rate maximizing interference alignment solution, solving ( [13] ) in its most general 
form is challenging. Consequently we pursue a numerical solution for certain special cases. 
Assuming K s = 1, (fT3]) simplifies to 



K a di j I , \ 

argmaxVVlog 1 + — . , (14) 

Fseo \ e n {CU + JWiHfcF.FraWj) e* n ) 

where subscript s refers to the single secondary user. As F s E O and R^ m (F s ) = K^ m (F s Q) for 
any unitary matrix Q, R^ m has the same value for all the points in the range of F s . Equivalently, 
we are searching for a d s -dimensional subspace of the vector space C M " and therefore, our search 



is confined to the Grassmann manifold. Hence, we can find locally optimum solutions for (14) 



using numerical optimization over the Grassmann manifold. We use the modified steepest descent 



in the complex Grassmann manifold algorithm (MGM) in [24| which requires the derivative of 



( [T4| ) with respect to F s and an initial guess. The problem in ( fT4| ) is not convex and MGM only 
guarantees convergence to a local optimum point. Therefore, a better initial solution will lead to 
a better final output of the MGM algorithm. Next we propose three initial solutions of varying 
degrees of complexity for the MGM algorithm. 
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A. Alternating Minimization Based Suboptimum Solution 

Let W = v^e n W 4 H !S , S = © © d i e n Q i e;, and d = K a d a . At high SDSfR, we can 



i=l n=l 



neglect the 1 and rewrite (14) as 



i=l n=l 



arg max — log det 
F s eo 



argmin log det ( I + -^S 1 W(I 



F s eo 



is 4 

7o d s 
To, 



Wfl; 



F S F*)W 4 



F s F*)W a 



(15) 
(16) 



where non- singularity of S (due to the rank constraint in (|2])) was used to derive ( fl~6] ) from ([15]). 
As lg ® F S F* = (Ij ® F fl ) (Ij <g> F*) and det (I + AB) = det (I + BA) for any two matrices 



A and B, (16) is equivalent to 



F s ee> 



argmin log det ( I + jilg ® F*)W*S^ i W(I 



(17) 



Let the eigenvalue decomposition of W*S 1 W be Uy^S^U*^. Note that rank ( W*S _1 W 



K a d a . If the rows of (Ij <g> F*) were equal to the linear combination of (M s — l)K a d a eigenvectors 
corresponding to zero eigenvalues in Sy^, the det in ( [T7] ) would attain its minimum value of 1. 
Therefore, ideally we would like to have 



I 



(18) 



where the columns of Uy^ are the columns of Uy^ corresponding to zero eigenvalues of 
W*S _1 W and A s is an (M s - 1) Ylf=\ <k x d i) d s combining matrix such that 



A„ A s 



A*TT*- TT - A 



(i,-®f:)(i 6 



Ij®F!F, = Ij®L 



I 



d s d' 



(19) 



Let U 



w 



.~.T 



r.T 



T 



, where ti q for q — , 

, Aj\, where Ap for p = 1 



1 d are row-blocks of dimension M s x 

d are column-blocks of dimension 



(M s -l)d. Let A s = [A, 
(M s — l)d x d s . Considering the block diagonal structure of Ij® F s , ( fj"8| ) can be rewritten as 

1 



U p A q = p,q 
UpAp = 1A p+ iAp + i p 



d, p^q 
-1 d-1 



(20) 



The system of equations in ( |20| ) can be further simplified to Usl s = 0, where a s = vec (A s ) and 
1A is found by vertically stacking the coefficients in the (M s d s )d 2 — (M s d s )d + (M s d s )(d — 1) = 

(M s d s ) ( d 2 — 1 ) equations given by pO] ). Note that as d 2 > M s , the nullspace of it only consists 



of the null vector and a non-zero a s exactly solving ( |20| ) cannot be found. Next we develop an 
alternating minimization algorithm to solve ( p~8| ) in the least squares sense. 

Assume F s in ( fT8| ) is given. We seek an A s to minimize ||Ij<8> F s — Uy^A s |||,. Formally, 



argmin tr 



h ® F. 



arg min tr ( 1^ 



f:f 5 +a:u;u^a s 



f:)u^a s 



(21) 
• (22) 
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Using ([T9]), d22l) simplifies to 



A, 



: arg min tr 



ha-** m 



f:)u w a. 



A s eo 



argmax 9te tr (A*U*^(Ij <g> F 



(23) 



where £He{.} selects the real part of a complex number. The optimization problem of ( |23| ) 
can be solved using the solution to the "Procrustes problem" p5| . Writing the singular value 
decomposition of U*^(I d ~ <g) F s ) as <fr 2 ] [Sf T ] T A*, the solution to ([23]) is given by 

A 



<&iA* with the maximum value of ( |2~3] ) equal to fHe (tr (Si)). 



Again consider (18) but assume that this time A s is given. We seek 



arg mm 

F s eo 



II.: 



U-wAs 



12 
\F 



arg max 9lt ( tr 



I 



)a:u^ 



Using the sparsity of Lj <g) F s , ( |24| ) can be written as 

F s = argmax 5He (tr (f s Efe=i A l U 



F a eo 



A; 



(24) 



(25) 



where «4.fc and U k are. defined in (|20j). Let * [U ± 0] [(A*) T (A;) T ] be the singular value 

decomposition of Ylt=i -^Mk- T ne solution to ( f25| ) is F s = Ai^* where the objective function 
attains its maximum value of D\t (tr (LTi)). Based on ([23]) and ( [25] ), the alternating minimization 
based approximate solution to ( [T8] ) is given by Algorithm [T] 



Algorithm 1 Alternating minimization based suboptimum solution 
Initialize F s 
repeat 

Find A s using ( [23] ) 

Find F s using ([25]) 
until Convergence 



After Algorithm 



converges, U-^A S 



d 

=i 



s and an initial solution for the MGM al- 
to the set of precoders in 



gorithm, approximately solving ( [14] ), can be found by restricting F 
{F l s , . . . , Fj 1 }. As the non-negative objective function of pT| ) reduces at each step of Algorithm 
[T] Algorithm [T] always converges. Convergence to a globally optimum point, however, is not 
guaranteed as O is not a convex set. Detailed results on convergence of alternating minimization 
algorithms similar to Algorithm [T] can be found in [26|. 



B. Interference Leakage Minimizing Suboptimum Solution 

The alternating minimization based initial solution suffers from high computational complexity. 
A less complex (and less accurate) initial solution to the MGM algorithm can be found by 
revisiting the constraint on the secondary user precoders outlined in ([6]). Instead of requiring 
the secondary transmitters' interference to be confined in the interference subspace of each 
active receiver, we minimize (in the least squares sense) the interference leakage caused by 

the secondary users through solving arg min ||H s F s |||n. A solution to this problem is given by 

F s eo 

setting the columns of F s to the d s least significant right singular vectors of H s . 
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C. DOF -Preserving Suboptimum Solution 



In both the initial solutions of Section V-A and Section V-B interference from the secondary 
transmitter is not confined to the interference subspaces of the active receivers. Thus, the DOF of 
the active users' network is zero. By aligning the secondary transmitter's interference at some of 
the active receivers, however, we can achieve a non-zero DOF and hence a better performance at 
high j . Specifically, as F s can not satisfy ([6]), we minimize the number of interference subspace 
dimensions at the active receivers to which the interference from the secondary transmitter is not 

aligned which translates to solving argmin ||H S F S || . In general, this is a combinatorial problem 

F s ee> 

(27j and exact solutions are out of the scope of this manuscript. Instead, we search for the largest 
subset of active receivers, K, a , such that the system of equations defined by PjH is F s = for 
all i e K, a has an exact solution (equivalently M s > ^2 di + d s ). In this fashion, F s will satisfy 

some of the linear equations in (|6]). Note that this solution is optimal when d a — 1. Moreover, 
although extending the alternating minimization solution to K s > 1 in not trivial, interference 
leakage minimization and DOF-preserving solutions can be adopted for arbitrary K s . 

VI. Numerical Results 

Consider an existing 3-user 2x2 MEVIO IA network where, through IA |9[ Section III. A], 
each transmitter sends a single stream to its corresponding receiver. For the rest of this section, 
we add nodes to this network and compare the total achievable sum rate of the resulting new 
network. All the channel and additive noise coefficients are distributed as CJ\f(0, 1). 

Single secondary user not interfering with the active users Assume K s = 1, M s = N s = 5, 
and d s = 1,2. The secondary user precoder can be designed using the optimum method presented 
in Lemma [2] or the random orthonormal and the selfish designs discussed in Section IV-A The 



average achievable sum rate of this network for various secondary user precoders versus 7 is 
shown in Fig. [2j For reference we show the sum rate when all the nodes cooperate through IA. 
As expected, the selfish precoding has the worst performance where, due to the uncoordinated 
interference, the DOF of the whole network is less than the original IA network. Moreover, both 
the optimum and random orthonormal precoders achieve 5 DOF, the maximum DOF attainable 
through IA (for DOF > 5 algorithm of [|9[ Section III. A] does not converge). At low j a , however, 
where interference and noise are not distinguishable, the selfish design has the best performance. 

Two secondary users not interfering with the active users Now, assume two M s x N s 
MIMO pairs (M S = N S ) are added to the network of active users. The secondary users can adopt 
selfish precoding, random orthonormal precoding, the self-optimizing precoding of Corollary 
[TJ or the iterative self-optimizing precoding of Corollary [2} The average total achievable sum 
rate versus M s for various precoders at two values of P a is shown in Fig. [3} The depicted 
upper-bound is for the case of no interference among the secondary users or between the active 

K a 

and secondary users. At 7 D = dB when N s — J2di is small, in contrary to the high j , the 

i=i 

selfish design is outperforming the other methods indicating the available spatial dimensions are 
better used to combat noise than interference. As discussed before, for M s = N s = K a + 1, after 
enforcing the zero-impact constraint on F& for k > K a , no further optimization of the precoders 
is possible explaining why self-optimizing, iterative self- optimizing and random orthonormal 
precodings have the same performance for M s = N s = 4. Random orthonormal precoding 
requires considerably less overhead compared to the other methods but achieves an acceptable 
performance. Also, the diminishing returns of increasing M s is more sever at low j a . 
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Three secondary users not interfering with the active users For the case of K s = 3 and 
M s = N s = 5, we consider selfish, random orthonormal, self-optimizing design of Corollary [T| 
iterative self-optimizing design of Corollary [2j and successive IA precodings. The average total 
achievable sum rate is shown in Fig. |4j In this case, as discussed in Conjecture [T] successive 
IA achieves the same DOF as IA over the whole network. The other methods, because of 
uncoordinated interference at the secondary receivers, cannot achieve more DOF than the original 
network of the active users. In other words, from Conjecture [TJ when M s < (K a d a + K s d s ) and 
d a = d s = 1, only successive IA can possibly achieve more DOF than the original IA network. 

Single secondary user interfering with the active users Consider the case where K s — 1 
and M s = N s = 3. Here, N s does not satisfy ([7]) and the new node is bound to degrade R^ m . 
Average R^m versus 7„ for the MGM algorithm initialized with the three solutions provided in 
Section [V] is depicted in Fig.[5j The alternating minimization algorithm approximately maximizes 
the sum rate of the active users and, as expected, the MGM algorithm converges to the best 
solution when initialized with it. In addition, the DOF preserving initial solution results in 
the best performance at high transmit power among all the initial solutions. Moreover, selfish 
precoding can be considered as a random initialization for the MGM algorithm and attains a 
good performance although the initial solution itself drastically decreases R^ m . 

Convergence of the MGM and the alternating minimization algorithms Both the MGM 
and the alternating minimization algorithms are guaranteed to converge to a (sub) optimum 
solution (see p4| Section VILA] and Section V-A). The convergence metric for the MGM 
algorithm is the magnitude of the descent direction on the Grassmann manifold and for the 
alternating minimization algorithms is the Frobenius norm of the change in F s at two consecutive 
iterations. The worst convergence rate of the MGM algorithm at two values of 7 Q initialized with 
the solutions of S ection |V-A"] and |V-B~| together with the slowest convergence rate of the alternating 
minimization algorithm over 1000 randomly generated channels are shown in Fig. [6] Note that 
the alternating minimization algorithm does not depend on j . Both algorithms converge in less 
than 500 iterations (the mean number of iterations of all the cases was less than 40). 



VII. Conclusions 

In this paper, we derived sum rate and DOF optimum/suboptimum precoders for a set of 
secondary users arriving to an already established network of active users performing MEMO 
IA. We showed how knowledge of the interference subspaces at the active receivers can be 
used to mitigate the performance degradation caused by the secondary users and how, in some 
networks, gains of IA are obtainable even if all the users do not perform IA together. We derived 
the minimum number of secondary transmit antennas, the zero-impact threshold, required for 
using the network resources without degrading the sum rate of the active nodes. 

When the secondary users have more than the zero-impact antennas, optimum and suboptimum 
precoders maximizing the achievable sum rate of the secondary users and the total DOF of the 
network were found for up to two and more than two secondary users, respectively. For up 
to two secondary users, we showed that a total achievable sum rate close to the upper-bound 
is obtainable. For more than two secondary users, the proposed successive IA techniques, for 
certain (and plausible) combination of number of secondary users and antennas, achieves the 
same DOF as if all the users had simultaneously performed IA. The successive IA method, as a 
decentralized IA technique, can be adopted in next generation wireless communication systems 
where different categories of nodes, such as base stations, femtocells and user terminals with 
different number of antennas and different levels of cooperation have to coexist. 
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When the secondary users have less than the zero-impact antennas, we proposed a precoder 
design based on steepest descent search on the Grassmann manifold. We derived three initial 
solutions for the numerical search, maximizing the sum rate of the active users at high SINR, 
minimizing the interference power leakage in a least squares sense, and maximizing the obtain- 
able DOF in the IA network. The sum rate maximizing initial solution was the most complex 
and resulted in the highest active users' sum rate. 

We assumed global CSI without considering the associated overhead. Also, with increasing 
number of antennas the interference mitigation power of the secondary nodes increases, but so 
do the overhead associated with training the additional links and relaying their CSI. Quantifying 
such overhead costs and the effects of a limited capacity feedback link is left for future work. 
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Fig. 1. Ka-usei MIMO interference channel where the ith transmit/receive pair is equipped with Mi transmit and Ni receive 
antennas. 



50 




Transmit Power (dB) 

Fig. 2. Average total achievable sum rate versus transmit power for varying precoder designs for the single 5x5 secondary 
node arriving to an already established 3-user 2x2 MIMO IA network. As expected, ignoring the active users results in the 
worst performance and reduces the attained multiplexing gain. Restricting the interference from the secondary transmitters at the 
active receivers increases the multiplexing gain of the network and can achieve the same DOF as if all the nodes had done IA 
together. Moreover, random orfhonormal precoding in the interference subspace of the active receivers achieves an acceptable 
performance close to the optimum design. 
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Fig. 3. Average total achievable sum rate versus N s at two secondary transmitters arriving to a 3-user 2x2 MIMO I A 
network for two SNR values. For small M„, the selfish design is optimum when the network is noise limited, but in the 
high SNR regime, restricting the interference at the active receivers considerably increases the sum rate. As expected, when 
M a = N s = K a + 1 = 4, all the interference avoiding designs have the same performance. Moreover, diminishing returns of 
increasing N s is more apparent at low SNR. 
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Fig. 4. Average total achievable sum rate versus transmit power for varying precoder designs at the three 5x5 secondary 
transmitters arriving to a 3-user 2x2 MIMO IA network. Ignoring the interference terms, although being optimal in the low 
SNR (and low interference-to-noise) regime |4|, reduces the total DOF of the network. Moreover, only successive interference 
alignment achieves the same DOF (6 in this case) as if all the nodes had jointly performed IA. 
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Fig. 5. Average total achievable sum rate versus transmit power for varying precoder designs for a single 3x3 secondary user 
arriving to a 3-user 2x2 MIMO IA network. As expected, MGM initialized with the alternating minimization solution achieves 
the best performance. The DOF preserving solution outperforms the other initial solutions at high transmit power while selfish 
and equal gain precodings result in the worst performance. 




Fig. 6. Worst convergence rate of the MGM and the alternating minimization algorithms over 1000 randomly generated 
channels. The convergence metric for the MGM algorithm is the magnitude of the descent direction on the Grassmann manifold 
and for the alternating minimization algorithm is the Frobenius norm of the difference between the F s at two consecutive 
iterations. Both the algorithms converge in less than 500 iterations in the worst case with the mean convergence rate of less 
than 40. 



